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This paper describes two attitude control laws suitable for atmospheric flight vehicles 
with a steady angular momentum bias in the vehicle yaw axis. This bias is assumed to be 
provided by an internal flywheel, and is introduced to enhance roll and pitch stiffness. The 
first control law is based on Lyapunov stability theory, and stability proofs are given. The 
second control law, which assumes that the angular momentum bias is large, is based on a 
classical PID control. It is shown that the large yaw-axis bias requires that the PI feedback 
component on the roll and pitch angle errors be cross-fed. Both control laws are applied to 
a vehicle simulation in the presence of disturbances for several values of yaw-axis angular 
momentum bias. It is seen that both control laws provide a significant improvement in 
attitude performance when the bias is sufficiently large, but the nonlinear control law is 
also able to provide improved performance for a small value of bias. This is important 
because the smaller bias corresponds to a smaller requirement for mass to be dedicated to 
the flywheel. 


I. Introduction 

This paper describes and demonstrates a Lyapunov-based nonlinear control law for an atmospheric flight 
vehicle that is “dual-spin” stabilized - that is, it includes an angular momentum bias along its yaw axis 
to gyrically stiffen its roll-pitch attitude dynamics. It has been observed * 1 that, if the bias momentum 
is sufficiently large, the vehicle’s roll-pitch dynamics approach those of an undamped second-order linear 
oscillator whose natural frequency scales with the magnitude of the bias. 

If this bias is sufficiently large - i.e. the system becomes gyrically dominant - the magnitude and 
bandwidth of control activity necessary to maintain a commanded attitude against disturbances can be 
dramatically reduced ; 1 furthermore, the approximately linear dynamics associated with gyric dominance are 
easily stabilized with a linear control law. 

Although a large yaw-axis angular momentum bias is desirable for pitch-roll stability and disturbance 
rejection, the fact that it is, in all probability, produced by a dedicated flywheel means that having a large 
bias tends to increase the vehicle’s mass. 

Designing a dual-spin atmospheric flight vehicle requires effectively trading the mass penalty associated 
with the flywheel against the advantages of being able to choose slower, less powerful fluidic control effectors 
for stabilization and command tracking. These advantages can be mundanely enhancing, e.g. cheaper or 
lighter control effectors, or they can potentially enable adoption of an emerging control effector technology - 
such as many current morphing technologies - that have difficulty in matching the authority or bandwidth 
available in traditional hinged control surfaces. 

The control issue that arises in selecting the “smallest possible” bias momentum flywheel for a given dual- 
spin-stabilized vehicle is that the gyric coupling in the dynamics is quadratic in its attitude rates, with this 
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nonlinearity only becoming insignificant if the bias momentum becomes dominatingly large; or alternatively, 
if it can be guaranteed that the vehicle’s attitude rates will be negligible, even in the face of disturbances. 
The consequence of this is that, for modest levels of angular momentum bias, a nonlinear flight control law 
is required. 

In the sequel, a Lyapunov-based nonlinear attitude control law for a rigid vehicle carrying a constant 
angular momentum bias is presented, and global stability proofs are provided for the case of time-varying 
commanded attitude, and for special cases in which the commanded attitude is constant, or in which there 
is feedback only of the attitude rates. These latter two simplifications provide insights that lead to a simple 
approach for selecting the gains. 

For comparison, a classical linear controller structure is described for the case where the vehicle is 
gyrically dominant (approximately linearizing the dynamics) with its angular momentum bias aligned with 
the yaw axis. Both controller structures are applied to the thrust-levitated research vehicle described in 
[2]. While undergoing substantial aerodynamic disturbance, the vehicle is controlled to execute a maneuver 
in which it rises vertically, then transitions to horizontal flight. The maneuver is simulated for identical 
disturbances with both linear and Lyapunov controllers, for several magnitudes of angular momentum bias. 
It is seen that the Lyapunov controller provides good performance with a very small angular momentum 
bias - corresponding to a flywheel system of no more than one percent of the vehicle mass - while the linear 
controller requires a dominatingly large angular momentum bias to provide comparable attitude performance 
- corresponding to a flywheel comprising up to five percent of the vehicle mass. 

The next section presents the Lyapunov-based controller and the vehicle dynamics for which it is intended. 
The stability theorems are given, along with a simple approach for selecting gains. Proofs of the theorems 
are given in Appendix A. Section 3 describes a classical controller for the gyrically dominant case. Section 
4 describes the particular dynamics of the simulated vehicle and disturbance environment for which the 
controllers are compared, and presents the results. Section 5 presents conclusions. 


II. Attitude Control via Lyapunov Approach 

The control laws in this study are intended for rigid body systems with dynamics that include a constant 
angular momentum bias, denoted by the symbol h B : 

v = -u> x v + C bo g 0 + — 4>f + v c - — Vs (1) 

m m 

J Co = — u; x (Juj + h B ) H- (j> r ~\~ t c — uj* h s — h s (2) 

where 

v c = — \c X UJ + c x u + u x c x uj] (3) 

m 

t c = c x C bo g 0 — c x v — c x v — Jut (4) 

Equations (1,2) are given in body coordinates, and v is translational velocity, m is mass, c is the vector 
from the body axes origin to the center of mass, w is vehicle angular velocity, J is the overall moment of 
inertia, and C bo g 0 is the gravitational acceleration vector transformed to body axes from NED coordinates. 
The control laws developed in the sequel assume that the angular momentum bias h B is assumed to be 
aligned with the 2 -body axis. The symbols (j>f and (f> T denote control forces and torques. In (3,4), r c and v c 
are secondary torques and forces, mostly associated with off-CG effects, where c is the vector from origin to 
center of mass. The quantities p s and h s are internal linear and angular momenta associated with moving 
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masses inside the vehicle. These generally would include momenta associated with rotating masses in the 
propulsion system, proof-mass dampers, and other such components. 

The first and most general Lyapunov controller in this section causes the vehicle to follow a commanded 
attitude history P(f) that corresponds to a angular rate command u>(i). These are related by the quaternion 
expression 


where 



P T P = 1 


[Pl^ 


D) 

-Pi 

-p 2 

-p 3 

V\ 

k'O 

-P 3 

^2 

Vi 

Fi 

T'O 

-Pi 

v 3 

-p 2 


VO 


= [p, B{v )\ , Ben 


4x3 


(5) 


( 6 ) 


Denoting the actual vehicle attitude as v, the tracking errors are 

Slo = uj — uj 
5u = v — v 

where 8v and Sui can be shown to satisfy 
where 
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Theorem 1 (Arbitrary Commanded Attitude) For the rotational equations of motion as given in (2), 
the feedback control torques 

4>t = 4>gyro + G u (uj — Uj) + G^fu — P) + Ju) — T c (11) 


guarantees that the tracking excursions relative to the commanded set of attitude and angular rates, v and 
uj, will be globally asymptotically stable if the gains satisfy the following conditions: 


(12) 

(13) 
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where 


4*gyro — ^ (JiJ H - ^ 5 ) 


(14) 


Noting that the inherent open-loop net gyroscopic torque on the vehicle is — w x ( Jco+hs ) as given in equation 
(2), the feedback torque, 4> gyr o, as given in equation (14) is a feedback linearization term. The proofs of 
Theorem 1 and subsequent Theorems 2 and 3 are outlined in the Appendix. Global asymptotic stablility 
means that Su> — > 03 X i, and 5v — > CUxi, as t — > 00 , for any initial condition dispersion about u and v. 

The next theorem restricts the control problem to tracking a fixed, rather than moving, commanded 
attitude. This permits the introduction of two tuning parameters: p and e in the 4> gyro expression from 
equation (14). For Theorem 2, 


<t>gyro{t,p) = £W X J LO + pU> X h B 


(15) 


where 


pe77:0<p<l (16) 

e G 77 : 0 < e < 1 (17) 


This permits the designer to separately scale the degree to which the control law is devoted to feedback 
linearization, via e, or to feeding back the gyroscopic torques associated with the angular momentum bias, 
via p. Note also that the now-constant P gives constant B which, in equation (12), results in constant G„. 


Theorem 2 (Fixed Commanded Attitude) For the rotational equations of motion as given in (2), and 


the scaled open-loop gyroscopic torque expression (15), the feedback control torques 

(J) T P) 4“ T Gj,(zz Zz) T c (1^) 

guarantee that the excursions relative to the commanded set of attitude, v, will be globally asymptotically 
stable if the gains satisfy the following conditions: 

G„ = /c„e77, > 0 (19) 

g u = [p-(i-p)A A ]j, Pen 3 * 3 , p< o (20) 

’ 0 1 0 1 

A\ = S\ 10 0 , 5a = i(A Pl - A P2 ) (21) 

_ ° ° ° J 

where (f> gyro > P , and e are given by equations (15, 16, 17), and X Pl and X P2 are 

Ap» — < = 1,2 (22) 

J i 


The X Pi of the above Theorem are interpreted by recalling 2 that yjX pl X P2 is the undamped natural frequency 
of the linear roll-pitch dynamics that occur when the system becomes gyrically dominant; that is, when 

\\M\/\\h B \\/^o. 
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As a control problem, Theorem 2 is a subset of Theorem 1. However, the restrictive set of control problems 
addressed in Theorem 2 allows a more flexible feedback control structure in terms of incorporating the scaling 
parameter for feedback-linearization e and the scaling factor p for feeding back the bias momentum induced 
gyroscopic torque. The control law in Theorem 2 guarantees that the vehicle attitude excursions relative to 
a constant commanded attitude, v, will be globally asymptotically stable. 


From Theorem 2, the net angular velocity feedback components in equation (18) can be expressed as follows: 


where 


<Mw) 


(frgyro T G U LU 

\P + €ui x — Q(p)] Jui 


Q(p) = p 


0 -A P2 0 
A Pl 0 0 

0 0 0 


+ (1 - p) A a 


We make the following observations: 


(23) 


(24) 


• The gain, Q(p), is a convex combination of a matrix resembling the state-space matrix of an undamped, 
roll-pitch-gyrically dominated system, and the sparse symmetric matrix, given by equation (21). 
If h B = 0, i.e., in the absence of bias angular momentum, Q{p) = 03 X 3- 


• If Ji = J- 2 , i.e., the vehicle is inertially axi-symmetric about its yaw axis, 


Ai — A 2 — A 0 


Aa — 0 3x3 


and 


/ 

Q(p) = p 

V 



(25) 


The design choice (p = 0, e = 0) will require no feedback torque proportional to both the gyroscopic 
torque due to the presence of bias angular momentum nor the inherent gyroscopic torque induced by 
vehicle angular velocities. This means that for a vehicle with bias angular momentum and which is 
inertially axi-symmetric about the yaw axis, a constant gain angular velocity feedback, 


cj) T {oj) = PJco , P < 0 (26) 

will ensure asymptotic stability (along with appropriate attitude error feedback). 

• The governing equation for the closed loop attitude response in Theorem 2 is 

JSu = —(1 — e)co x Juj — (1 — p)uj x h B + G u Suj + G v 8v (27) 

The ubiquitous kinematic nonlinearity is evident. On the other hand, only the value of the design 
parameter, e, affects the presence of nonlinear dynamics in the closed loop response. In particular, 
only the value of e = 1 ensures that there is an exact feedback linearization of the inherent nonlinear 
gyroscopic torque, i.e., the closed-loop system responds dynamically nonlinearly if e 7^ 1. However, 
from the control law given by equation (18), notice that only e = 0 will ensure a linear angular velocity 
feedback. Hence for the range 

0 < e < 1 (28) 

both the control law and the corresponding closed-loop response involve nonlinear dynamics. A conse- 
quence of this is that linear stability analysis and control design methods may yield unreliable results. 

• This special form of the control law is analogous to previous results such as [3] , [4] , [5] , and [6] . 
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Consider stabilizing the vehicle attitude using only angular velocity feedback. Obviously without attitude 
feedback, the vehicle’s attitude cannot be tracked to any commanded value but it can be regulated to some 
constant. Theorem 2 applies with the added condition of zero attitude feedback gains, G v = 03 X 4 - However, 
notice that the negative definite condition on the key angular velocity feedback gain parameter, P, does 
not provide a clear guidance for designing this gain matrix for a physical application. For this reason, we 
consider an additional condition on this gain matrix, namely, we limit the gain parameter, P, to a diagonal 
structure. The mathematical structure of the resulting feedback- linearized closed loop system is analogous 
to the damped open-loop system as studied in some detail in [ 2 ] for sizing the requirement on angular 
momentum bias. This analogy provides added insight for controller design. 


Theorem 3 (Angular Velocity Feedback Only, P = diag(ai,a2,a 3 ) < 0 ) For the rotational equations 
of motion as given in ( 2 ), the feedback control torque 

(f> T = figyro^p) T G^OJ T c (29) 


guarantee that the vehicle dynamics in ( 2 ) will be globally asymptotically stable (u> — > O3X1 as t — > 00) if the 
feedback gain matrix satisfies the following: 


where 
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(30) 

(31) 

(32) 

(33) 


(34) 

(35) 


The 4> g yro{p ) term in above control law fully feedback- linearizes the dynamics, resulting in the linear, time- 
invariant closed loop system: 

w = J~ 1 AJu (36) 

where 


A 

7 

A 


ai —7 0 

7 o 2 0 
0 0 03 


(l-p)A 

^(-Vi + 2 ) 


(37) 

(38) 

(39) 


Since eigenvalues are invariant with respect to a similarity transformation, the eigenvalues of the closed loop 
system in equation (36) are: 


Ai ,2 = -< ± joj n x/l - C 2 (40) 

A 3 = a 3 (41) 
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where £ and uj n denote the damping ratio and natural frequency of the roll-pitch closed loop dynamics and 
are related to the design parameters by 

= -^(ai + 02 ) (42) 

<4 = 7 2 (43) 

As stated above, Theorems 2 and 3 are simplified restrictions of Theorem 1 that provide additional insight 
for control design: Theorem 2 for attitude tracking and trades on feedback linearization and gyroscopic torque 
compensation, and Theorm 3 for attitude stabilization. Focussing on the latter issue, note from equations 
(39, 42, 43), that Theorem 3 gives three parameters that can be manipulated to adjust the natural frequency 
and damping ratio of the roll-pitch dynamics: a±, a 2 , and p. These parameters are selected by making two 
assumptions. First, assume that only the closed loop frequency response is of interest, rather than mode 
shapes. Secondly, invoke physical intuition to assume that the least amount of control effort will be expended 
if the control law is dedicated to providing damping, while leaving the closed loop frequency equal to its 
open loop value. The first assumption and equation (42) give 


Quin = —a ai = a ,2 = a 
The second assumption and equation (22) give 


— 


7M4 


(44) 


(45) 


and from equation (43), 


P= 1- 


2 \J J\Ji 
J\ + J 2 


(46) 


In summary, to obtain desirable closed loop transient properties, p can be chosen from equation (46) to 
attain a natural frequency of u n while ( 01 , 02 ) can be chosen to satisfy equation (42) for a selected damping 
ratio, (. The simplest choice for these latter is to set them equal, per equation (44). 


III. Attitude Control via Classical Feedback 


In order to obtain simple heuristics for control the following reasonable assumptions are made 


c ft 0 1 

\\M\ « ||M IIM « IIM \ 

\\M\ « \\M\ IIM « \\M\ ) 

With the above assumptions, the following approximate equations are obtained: 

v « -u] x v + C b0 g 0 + — (j) f 

m 

J uj ft — co x (J to T /i B ) T (j)r 


(47) 


(48) 

(49) 


These approximate equations are used only for the purpose of deriving control laws; all closed loop simulations 
shown in this study are based on the full equations. 
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Consider an attitude control law based on classical feedback ideas. The approximate attitude dynamical 
equations, for an air vehicle which has a significant bias angular momentum, i.e., gyrically dominated in the 
body z-axis, as given by equations (49), are rewritten in the following expanded form: 

(50) 

(51) 

(52) 


LOi « 

— \ pi id2 + — 
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d>2 ~ 

\ , T 2 
A P2 wi + — 
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where 


- (ti,T 2 ,T 3 ) 


(53) 


A. Roll-Pitch damping by angular velocity feedback 

Consider the feedback control law: 

j- = k Dii uj l + k Di2 u 2 + u 1 
Y = k D2i ui + k D22 UJ2+U2 


(54) 

(55) 


where (ui, U 2 ) are attitude error dependent feedback terms used subsequently to track commanded attitudes. 
To determine the gains, (k D ,k D ,k D ,k D ), consider the closed-loop dynamical equations: 
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The above closed-loop roll-pitch dynamics clearly correspond to the damped open-loop dynamics as given 
by equation (3.1) found in page 7 of [2]. In particular, (k Dii ,k Dnn ) correspond to damping coefficients, 
(k D ,k D ) correspond to stiffness changes modifying the precession frequency, and ( 111 , 112 ) correspond to 
external torque disturbances. To mitigate control requirements, we do not wish to change the precession 
frequency of the gyrating open-loop system so that we choose 


= 0 


(57) 


Additionally, we ignore the “external torque disturbances” so that the roll-pitch closed-loop dynamical 
equations given in equation (56) reduces to the following: 
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The closed-loop eigenvalues of the above system are given by (see equations (3.5) to (3.8) on page 7 of [2)] 
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At this point we can select a desired damping ratio and then solve for the gains, k D and k Doo , from equations 
(60) and (61). However, as a simpler baseline approach, consider the following additional assumption: 

= k D22 = k (63) 

With this simplifying assumption, equations (60,61) give the following design condition: 

k = —\ 0 —^== < 0 (64) 

In principle, one can choose an arbitrary closed-loop damping ratio without changing the precession frequency 
of the open-loop system by choosing the angular velocity feedback gains that satisfies equations (60, 61, 57). 
Based on equation (64), a design rule of thumb could be inferred for the following nominal condition: 

- -f ~ C; C < 0.3 (65) 

which indicates that the ratio of the angular rate gains to the precession frequency is indicative of the vehicle 
damping ratio. 


B. Attitude tracking by attitude error feedback 


With the asymptotic convergence of uq and w 2 guaranteed by angular velocity feedback as described earlier, 
the requirement for convergence to commanded attitude angles still remains. This tracking is achieved by 
feeding back attitude errors, as outlined next. 


Suppose 


t 2 = 0 


so that the pitch-axis dynamics from equation (51) become 


W 2 ~ A p2 0q 


( 66 ) 

(67) 


By taking the time derivative of pitch equation (67), substituting (50) for cLq, and using the roll axis control 
law given by equation (54), and substituting for uq using equation (67), it can be shown that 

w 2 - -r^-k Dii u) 2 + (A P2 A Pl - \ P2 k Di2 )u> 2 ~ A P2 ui (68) 

We make the following observations on equation (68): 


• The pitch angular velocity state is uncoupled when viewed as a second-order oscillator. 

• The pitch angular velocity is driven by roll torque. 

• The coefficients of the second-order oscillator are such that any transient free response will be asymp- 
totically stable. 


Since the asymptotic convergence, w 2 — > 0 as t — > 00 , does not ensure convergence to a commanded pitch 
angle, the additional control input torque component, iq, is introduced. So assuming the small angle linear 
kinematics for pitch axis rotation, 

O 2 = w 2 (69) 

it can be shown that a PI feedback of the pitch angle error 


u 1 = - 




( 70 ) 
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will guarantee convergence to any step command in pitch angle, 02 md . Notice that the pitch angle error is fed 
back as a commanded torque in the roll axis, which is intuitively consistent for a gyrically dominant system 
with the bias momentum in the yaw axis. However, design details for choosing gains for the PI control law 
are not addressed in this study given the availability of various techniques based on classical control design 
ideas (see for example Chapter 3 in [7]). 

Using the above approach but assuming t\ = 0, the required PI feedback to track commanded roll angle 
takes the form 



The roll angle error is fed back as a commanded torque in the pitch axis. As for tracking yaw angles, the 
bias momentum which strongly cross couples roll and pitch axes do not affect the yaw dynamics except 
through the nonlinear vehicle gyroscopic component ^ ' ,2 J dl ^ W 1 W 2 , which is a small, compared to the wheel 
gyroscopic torques, X Pl L 02 and A P2 0 Ji- Hence tracking the yaw angles will involve a direct application of PI 
feedback of the yaw angle errors. 


To summarize, the feedback control law for providing damping and attitude tracking is as follows: 

k. 
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(72) 

(73) 

(74) 


The yaw axis control law is composed of a feedback linearization term, a linear feedback of yaw angular 
velocity to provide damping and a PI feedback of the yaw angle error for tracking. The roll and pitch control 
laws consist of two components; a linear feedback of roll and pitch angular velocities to provide damping and 
PI feedback of the roll and pitch angle errors for tracking. However, the tracking angular errors are cross- 
coupled. The gyrically dominant influence of the bias angular momentum placed along the vehicle body yaw 
axis essentially feedback-linearizes the roll-pitch dynamics. Hence there is no need to use nonlinear feedback 
control for linearization. Finally, note that damping ratios on the precession frequency can be arbitrarily 
placed using for example, equation (64), while not changing the open loop precession frequency. 


C. Uncoupled small angle attitude kinematics 

With this single-loop approach, the angular tracking error terms are treated independently and conveniently 
so that no rotation sequence need to be considered. On the other hand, the kinematic nonlinearity inherent 
for large angular motion (say beyond 30 degrees deviation) limits the applicability of this control law. 
However, notice that the angular rate and acceleration feedback terms still apply for large angular motion so 
that a judicious choice of control gains may result in satisfactory closed loop performance. For small angular 
motions, these attitude tracking errors can be approximated directly from the Euler parameters as follows 
(see for example Chapter 2 in [8]): 


02 

= arcsin [—Cbo i3 (v)} « 2v 2 

(75) 

03 

= arcsin[Cbo i2 (is) / cos 0 2 ] ss 2v\ 

(76) 

01 

= arcsin [Cbo 23 ( v)/ cos 0 2 ] « 2v 3 

(77) 


Notice that the above equations relate 321-sequence Euler angles (in radians) to the unique direction cosine 
matrix, which is parameterized by a set of Euler parameters. As is well known, the above relationship suffers 
a kinematic singularity when the pitch angle equals ±7r/2. 
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IV. Control Law Comparison via Simulation 


This section describes numerical simulation experiments conducted to characterize the performance of 
the linear and nonlinear (Lyapunov) controllers for a range of magnitudes of the angular momentum bias 
h B . As mentioned in the Introduction, it is desirable to employ the smallest h B that will adequately stiffen 
the vehicle’s attitude dynamics, since that results in the smallest flywheel mass. In this study, the tradeoff 
between vehicle mass and closed-loop performance is realized by posing several flywheels with differing 
mass but common physical assumptions: All flywheels were assumed to be thin steel rings with tensile 
strength of 80kpsi and a 30krpm rotational rate, with their radii chosen so that the hoop stress would be 
80kpsi. Four cases were considered, indexed on the flywheel’s fraction of the total vehicle mass. These were 
m.fiywheei = {0, 1%, 5%, 20%} • m ve hide- For each of these four cases, a classical controller and a nonlinear 
controller were formulated, and the vehicle was controlled to follow a commanded trajectory in the face of 
aerodynamic disturbances. 

Figure 1 shows a sketch of the simulated vehicle used for the control comparison experiments, along with 
its key nominal parameters. It is a small square platform with a shrouded propeller at each corner. The 
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Parameter 

Symbol 

Value 

Units 

Vehicle mass 

m 

20 

kg 

Vehicle length, width, height 

h, h, h 

.9, .9, .3 

m 

Vehicle second moment of inertia 

h, h, h 

.59, .58, 1.15 

kgm 2 

Ducted fan rotor inertia (each) 

ia F 

.002 

kgm 2 

Ducted fan speed (range) 

ids 

5 to 12 

krpm 

Ducted fan control vane (range) 

Q v 

—30 to 30 

deg 


Figure 1. Schematic of a generic quad ducted fan with bias momentum wheel. 


vehicle’s body axes have their origin at the center of mass, and are oriented as shown more clearly in Figure 
2, with the 2 -body axis pointed into the page. Thrust vectoring vanes are mounted downstream of each 
propeller, and each set of vanes diverts the thrust along the body axis normal to the body axis that passes 
through the thruster. Figure 2 displays both the vehicle body axes and four local coordinate systems for the 
four thrusters: Fp . , j = 1,2, 3,4. The thrust deflections due to thrust vectoring vanes are parallel to the 
y axes of the Fp. frames. The momentum wheel unit is assumed to be rigidly attached to the center of the 
vehicle which, itself, is assumed rigid. The spin axis is assumed to be collinear with the vehicle yaw axis. 

The commanded trajectory to be followed by the vehicle is shown in Figure 3, and consists of a vertical 
ascent that transitions to level cruise. Starting from a hovering position, the vehicle is commanded to climb 
at a flight path angle of 80 degrees while accelerating for a duration of 10 seconds. This is followed by a 
transition of flight path angle from 80 degrees to 0 degrees in 5 seconds while continuing to accelerate. At 
the end of 15 seconds the vehicle is commanded to a level flight while the airspeed continues to accelerate 
for the next 5 seconds. From 20 seconds into the flight, steady level cruise is commanded for the next 20 
seconds. Note that for most of this trajectory, the commanded attitude is constant. 
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Figure 2. Top view of ducted fan frame orientations relative to platform body frame. 


The air disturbances to which the vehicle is subjected during flight were assumed to come from three 
sources: (1) basic aerodynamic loads as driven by the steady wind states (V T ,a,(3,u>), (2) air turbulence as 
defined by a Dryden model as described for example in Chapter 4.5 of [9], and (3) ground effects. The basic 
aerodynamic loads on the vehicle are shown in the left side of Figure 4. Notice the steep increase in drag 
forces and the pitch down aerodynamic torque as the vehicle increases to cruise speeds. The net aerodynamic 
loads with the additional simulated turbulence as described by a Dryden model with the following turbulence 
parameters: intensity a w = 3.8 m/s, scales (L x ,L y , L z ) = (121, 60, 30)m, are shown in the right side of Figure 
4. To simulate the air loads on a vehicle that takes off vertically from a ground position and then cruising at 
low altitudes which may be gusty due to ground obstacles such as buildings and uneven landscape, random 
gust torques having a bandwidth of 3 Hertz in all axes are assumed. The amplitude of this random signal is 




Time (sec) 


Figure 3. Commanded vehicle trajectory: airspeed (m/s), path angle (deg), heading (deg). 
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no air disturbance air disturbance 




no air disturbance 


air disturbance 




Figure 4. 


Aerodynamic Loads on the Vehicle About Trim: without turbulence (left), with turbulence (right). 
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set to taper off linearly from ground to an altitude of 5 meters where it is set to zero above it. In addition, 
low level random disturbance torques are added to account for persistent gusts in low altitude atmosphere. 
Figure 5 shows the simulated random gust air loads due to the ground effects and low altitude flight. 


air disturbance 



Figure 5. Random gust loads on the vehicle due to ground effects and low altitude flight. 


Figures 6 and 7 summarize the closed-loop performance of the vehicle in following the control trajectory 
in terms of rms values; for Figure 6, the deviation of attitude and velocities from commanded values, and for 
Figure 7, deviations of control variables from their trim settings. The rms values are computed as the square 
root of the time integral of the squared deviation, divided by the integration interval. The “case numbers” 
labelled on the horizontal axes of the plots refer to the combination of control law type and flywheel (bias 
momentum) size. The prefix “C” refers to a classical control design and “L” to a nonlinear one. The numbers 
0, 1, 2, 3, refer to 0%, 1%, 5%, and 20% flywheel mass fraction, respectively. 

The first thing to note in Figure 6 is that all the bars for case Cl extend to the top of the charts. This is 
because the authors were unsuccessful in finding a stabilizing linear controller with the structure described 
in Section 3. The 1% flywheel for this case did not generate enough h B to gyrically dominate the dynamics, 
but did generate enough to invalidate the linear model assumed by the controller. The case CO, with no 
bias momentum, actually worked better. In all cases, the “L” controllers did somewhat better in terms of 
rms behavior than the corresponding “C”s, particularly in the case of LI versus Cl. The second thing to 
note in Figure 6 is that the performance of the linear and nonlinear control laws is quite comparable for the 
larger flywheel masses, though the Lyapunov controller was marginally better. In Figure 7 note that the 
Lyapunov controller appears to be somewhat more aggressive than the classical controller in its use of fan 
speed, though the controllers are otherwise comparable, at least at the 5% and 20% flywheel mass fractions. 
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o of Attitude Error (NED) 



o of Angular Velocity 



a of Ground Velocity (NED) 



Case Number 


Figure 6. vehicle response comparisons. 
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a of Vane Deflection a of Vane Deflection Rate 




a of Fan Speed o of Fan Acceleration 




o of dfan Forces 



Case Number 


a of dfan Torques 



Figure 7. control effort comparisons. 
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V. Conclusions 


This paper has described two controller structures suitable for an atmospheric flight vehicle that includes 
a steady yaw-axis angular momentum bias. A nonlinear controller based on Lyapunov stability theory was 
described and its stability was proven, for the full formulation, and for two restricted cases that allowed a 
simplified gain selection approach. In addition to this, a classical PID-type linear controller structure was 
described for the case where the vehicle’s attitude dynamics was gyrically dominant - where the angular 
momentum bias was sufficiently large that the nonlinearities in the dynamics were masked by the bias. It 
was shown that this controller structure requires that the PI roll and pitch error feedback must be cross- 
fed. Both controllers were demonstrated by applying them to a simulation of a thrust-levitated vehicle 
that is commanded through a maneuver in which vertical flight transitions to horizontal in the presence of 
aero disturbances. The controllers were simulated for several levels of angular momentum bias, and it was 
seen that, while both controllers could provide good response for substantial angular momentum biasses, 
the nonlinear controller was necessary for controlling the system for a very small angular momentum bias. 
This was important because, in the quest for the lightest possible vehicle, there is considerable pressure to 
minimize the mass of the bias momentum flywheel. 


A. Proof of Theorems 


The proof is based on Lyapunov stability theory (for more details see for example [10]) and shows the 
asymptotic stability of the equilibrium point, (Su> = 03 X i, 8u = 04 X i). The theorem assumes the dynamical 
system given by equations (2) whose attitude error kinematics are described by equation (9). The closed 
loop dynamical equations can be obtained by substituting the feedback control torque in equation (11) into 
the dynamical equations (2). However, to simultaneously prove Theorems 1 and 2, let us assume the more 
general form of <j) gyr o as given in equation (15). This results in the closed loop response equations 

JSuj = — (1 — e)u> x J (jj — (1 — p)io x h B + GujSuj + G v 8v (78) 

where G„ is given by equation (12) and G u by equation (13). 


For the above closed loop system, consider the following function about the null equilibrium point: 




*3x3 


k v L 


i/*4x4 


Sh 

Sv 


where 


5h 4 J5uj 


The above choice of function satisfies the necessary Lyapunov function properties: 


P(0,0) = 0, and P>0,V 


Sh 

8v 


± 0 


(79) 

(80) 

(81) 


Notice that the variable 5h in the previous equations can be substituted by Slu without any change in the 
property of V since J is nonsingular. 

It can be shown that along any solution trajectory of the system defined by equations (78) and (9), 

V = Sh T [— (1 — e)w x Jlo — (1 — p)ui x h B + G u Su> + G v 5v\ 

( 0 


1 T 

+—k v Su T 


Suj 


M Sv 


( 82 ) 
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By partitioning and defining as given in equation (9) we can rewrite 


Suj 


= BSco 


In addition, since [w] is skew-symmetric 


Sv T [u>\Sv = 0 


and furthermore, since the attitude feedback gains satisfy equation (12), the terms 

5h T G„Sv + ^k v 5v T BSui = 0 


(83) 

(84) 

(85) 


Hence, equation (82) 21) and A is defined by equation (39), to rewrite the V expression in equation (??) as 
follows: 


V = -Sh 1 



f Ml 

(1 - e)u - (1 ~ p) 

0 


l a )\ 


Jlu + (1 — p)8h 1 A\ Juj + Sh 1 G u J L Sh 


( 86 ) 


At this point, the proof for Theorems 1 and 2 are considered separately. 


Theorem 1: e = 1 and p = 1, i.e., the complete feedback linearization of all gyroscopic torques is assumed. 


The V expression in equation (86) reduces to 

V = 5h T G ul J~ 1 Sh (87) 

To satisfy Lyapunov stability condition we note that in addition to conditions in equation (81): 

V = 6h T G U) J~ l Sh < 0, W(Sh, Sv) (88) 

<*= V < 0, WSh^O (89) 

o G u J~ l < 0 (90) 

& 3 G u : G u = PJ, Pe U 3x3 , P < 0 (91) 

The condition for the gain matrix, G^, in equation (91) is identical to that of equation (13) found in 
Theorem 1. This proves the Lyapunov stability implicit in Theorem 1. To show that the system is actually 
asymptotically stable, note from equation (88) that 

V = 5h r G u J- 1 5h = 0 iff Sh = 0 (92) 


This suggests a proof of asymptotic stability by showing that the only feasible solution in the subset of 
solutions (Sh = 0, 5v) satisfying V = 0 is the null solution, which is our equilibrium solution of interest. By 
substituting these additional conditions into the closed loop system equation (78) and noting the requirement 
k v > 0, it follows that 

B T Sv = 0 3x i (93) 

It follows from equation (9) that a basis vector for the null space of B T in the above equation is v £ lZ 4xl . 
This means that any solution, Sv to the above equation is a scaling of v itself, and since an Euler parameter 
vector is unit norm bounded, we can conclude that Sv = CLxi- Hence the only solution that satisfies 
V = 0 condition is the null solution (Sto = Oaxi,^ = CLxi). By a corollary of the Lyapunov theorem 
(see for example, [10]), this proves asymptotic stability. Finally, observe that with the above choice of 
the Lyapunov function which is continuously differentiable, radially unbounded, and positive definite, the 
asymptotic stability result applies globally as well. 
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Theorem 2: to = 0 but for any values of e € [0 : 1] and p £ [0:1]. This implies that 


Sh = JSuj = Jlo 


(94) 


which represents the set of control problems involving attitude regulation. Using equation 94, the V expres- 
sion in equation 86 reduces to 

V = 5h T [(1 — p)A\ + G u J~ l \ Sh (95) 


since 


Sh 1 


To satisfy Lyapunov stability condition note that: 



(Ml 

(1 - e)w - (1 - p) 

0 


(Ml 


Sh = 0 


(96) 


V = Sh T [(1- p)A x + G u J~ 1 ]Sh<0, V(/i, Si/) (97) 

<£= U<0, VSh^O (98) 

(1 - p) Aa + G^J- 1 < 0 (99) 

3Ga, : = [P — (1 — p)A\]J, P e U 3x3 , P< 0 (100) 


The condition for the gain matrix, G u , in equation (100) is identical to that of equation (20) in Theorem 
2. This proves the Lyapunov stability implicit in the Theorem. To show that the system is actually 
asymptotically stable, note from equation (97) that 

V = 6h T [(1 - p)A x + G^J- 1 ] Sh = 0 iff Sh = 0 (101) 

As in the asymptotic stability proof in Theorem 1, this suggests a proof based on showing that the only 
feasible solution in the subset of solutions ( 6h = 0, 5v) satisfying V = 0 is the null solution. Hence the 
remaining steps involve substituting the particular additional conditions in Theorem 2 into the closed loop 
system equation (78). The end result is equation (93) with the conclusion that the only solution that satisfies 
V = 0 condition is the null solution. The reasons for global asymptotic stability holds here as well. 


Theorem 3: Theorem 3 is a special case of Theorem 2. 
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